Mixed Ab-initio Quantum Mechanical and Monte Carlo Calculations 
of Secondary Emission from SiC>2 Nanoclusters. 
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A mixed Quantum Mechanical/Monte Carlo method (QMMC) for calculating Auger spectra from 
nanoclusters is presented. The approach, based on a cluster method, consists of two steps. Ab-initio 
quantum mechanical calculations are first performed to obtain accurate energy and probability 
distributions of the generated Auger electrons. In a second step, using the calculated line-shape 
as electron source, the Monte Carlo method is used to simulate the effect of inelastic losses on 
the original Auger line-shape. The resulting spectrum can be directly compared to 'as-acquired' 
experimental spectra, thus avoiding background subtraction or deconvolution procedures. As a 
case-study, the O K-LL spectrum from solid Si02 is considered. Spectra computed before or after 
the electron has travelled through the solid, i.e. unaffected or affected by extrinsic energy losses, 
are compared to the pertinent experimental spectra measured within our group. Both transition 
energies and relative intensities are well reproduced. To the best of our knowledge, this is the first 
time that Auger spectra are computed ab-initio, with no free parameters involved. 
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I. INTRODUCTION 

In recent years, remarkable developments in two 
branches of semiconductor physics have been achieved. 
One is the theoretical advance in the study of band gaps 
and quasi-particle spectra in systems, notably transition- 
metal oxides, where electronic correlation creates states 
other than simple Fermi gas [l[ . The other is the science 
of nanomaterials, whose electronic, spectroscopic and 
transport properties bridge the gap between molecules 
and condensed matter (2j. The search for new technolo- 
gies in nanoelectronics based on silicon and carbon nano- 
materials [H [i[ has been fueled by the potential applica- 
tions of such nanoclusters as junctions in optical telecom- 
munication networks or in integrated circuits where se- 
rious obstacles are represented by both the signal speed 
and spatial density. Despite many recent advances in the 
development of carbon based electronics [5( , silicon still 
represents the most widely used physical means. Further- 
more, a miniaturization of silicon based devices would 
retain all our understanding of the physics of this semi- 
conductor and the technological background of the mi- 
croelectronics industry. But while the knowledge of the 
electronic properties of silicon and silicon oxides has ad- 
vanced greatly, largely due to the development of theories 
capable of treating the correlation energy Q , theoretical 
results on the secondary emission and energy loss mech- 
anisms are still the subject of discussion [7[. 

In this paper, we focus on the O K-LL spectrum from 
SiC>2 clusters. We present a combination of an ab-initio 
Quantum Mechanical method with Monte Carlo simu- 
lations (QMMC) to calculate non-radiative decay spec- 
tra from Si02 nanoclusters of different size taking into 
account post collisional interactions (PCI) due to elec- 
tron dynamical screening. The method is validated by 
comparison to measured spectra obtained in our group. 



In this regard, we move along two independent direc- 
tions. On the one hand, we start from the measured 
spectrum, affected by all kinds of energy losses suffered 
by the electron on its way out of the solid, and we recover 
the 'true' Auger spectrum by conventional deconvolution 
procedures. On the other hand, using the computed O K- 
LL spectrum as electron source, we use the Monte Carlo 
method to simulate the effect of inelastic scattering on 
the original electron energy distribution. Key ingredi- 
ent for both procedures is the inelastic scattering cross 
section or, equivalently, the Reflection Electron Energy 
Loss Spectrum (REELS). While the well studied @,[l,[iQl 
O K-LL spectrum from Si02 clusters is taken here as a 
case study, we believe that our approach is sufficiently 
general to be applied to the interpretation of core elec- 
tron spectroscopies, such as Photoemission (PES), Auger 
(AES) and Autoionization by taking into account dy- 
namical screening effects in a variety of nanoclusters and 
condensed matter systems. 

In order to gain new insights into the nature of correla- 
tion and of continuum states in nanoclusters, a number of 
well-defined questions have to be addressed. What is the 
nature of the fundamental process occurring in a scatter- 
ing experiment, where an electron or a photon is injected 
in the topmost layer? What is the nature of the core- 
hole and double-hole states when electronic relaxations 
and charge disorder are included? How does hole local- 
ization distort the double-hole density of states (DOS) 
and manifest itself in the scattering process? What is the 
role of dynamical screening in the Auger process? Can 
the above considerations be combined to understand ex- 
perimental data in a unified framework? Are results size 
dependent? 

Many models, described in several reviews [HI [H, 
EH 13, have been suggested to tackle the calculation 
of Auger spectra in solids, a difficult task due to the 
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high number of energy levels involved and the contri- 
bution of additional degrees of freedom which broaden 
lineshapes. In particular, computational and theoreti- 
cal procedures are represented by methods developed by 
Cederbaum 0hn0 E3 and Cini-Verdozzi []J|. In 

each of these cases, Auger decay probabilities were com- 
puted by means of atomic Auger matrix elements with 
the many-body Green function projected onto atomic 
states and the spectral lineshape obtained by calculating 
the residua of such a projected quantity. This approach 
revealed very valuable in the calculation of atomic and 
molecular [13, [2(| HH Auger spectra but the extension 
to solids is fairly difficult and poor in results. Verdozzi 
and Cini [TH, [22j improved such an approach to sepa- 
rate the atomic features from those characteristic of a 
periodic structure (see below). The striking difference 
between our approach and those used by all the above 
cited authors (and further extensions of their work) lies 
in the computation of the many-body Green function. In 
fact, in our method the many-body Green function is pro- 
jected onto 'localized' states, which are a 'local mixture' 
of atomic states. Therefore, differently from previous 
works, our local projectors are multicentered and may 
include some atomic space points. In other words, the 
role played by our projectors is somehow similar to that 
played by localized Wannier functions in solid state calcu- 
lations. It is clear that the use of 'localized' (not atomic) 
Auger matrix elements will prove to be very valuable 
in the treatment of transitions, such as valence-valence- 
valence (V-VV), where initial states are not atomic, but 
delocalized in nature. The quest of multicenter projec- 
tors in the calculation of Auger matrix elements is of 
paramount importance, because the use of atomic Auger 
matrix elements, acceptable in molecular spectra, would 
result in a crude approximation for condensed matter 
core-hole spectroscopy. 

As for core-valence-valence (C-VV) spectra, their un- 
derstanding is presently based on the use of atomic tran- 
sition matrix elements and it is rooted in Lander's orig- 
inal idea [23| that such spectra should basically reflect 
the self-convolution of the valence band DOS, so called 
'band-like' spectra. While this is true for solids such as 
Li and Si, it fails for transition metals, Zn for example, 
whose C-VV spectra are instead referred to as 'atomic- 
like'. Recognizing that correlation effects between the 
two final-state holes were responsible for such spectra, 
led to the Cini-Sawatzky theory [24], [Hj] which, derived 
at first only for systems with initially filled bands, has 
since become the framework for any quantitative inter- 
pretation of C-VV spectra. Within this model, hole-hole 
correlation energy is a semi-empirical parameter to be de- 
termined by comparison with measured spectra. In par- 
ticular, the effect of hole-hole repulsion on the O K-LL 
spectrum from bulk Si02 has been recently investigated 
by Van Riessen at al. Busing the Cini-Sawatzky model 
extended by Ramaker [9(. The spectrum is computed 
by multiplying experimental intra-atomic transition ma- 
trix elements times the distorted DOS due to hole local- 



ization. Despite the O K-LL spectrum could be satis- 
factorily reproduced, matrix elements and characteristic 
parameters were taken from atoms, a fact that can be 
inappropriate for ionic crystals, such as Si02, where oc- 
cupation numbers are different from the atomic ones. 

To the best of our knowledge, the calculation presented 
here represents the first attempt to address the ques- 
tions listed above from first principles, with a full ab- 
initio treatment of hole-hole correlation and taking into 
account collective excitations without free parameters, 
except for a constant factor to normalize spectral inten- 
sity. We will focus on the so-called principal or normal 
Auger process, whereas the contribution of satellites (as- 
sociated with shake processes) to the Auger spectrum is 
neglected. 



II. THEORY 

The general theoretical framework for interpreting 
Auger decay in nanoclusters is the theory of scattering 
[26| . where the initial state consists of a projectile, in 
our case typically photons or electrons, and a target in 
its ground state, while the final states are characterized 
by an ionized system and one or more electrons asymp- 
totically non interacting. The main problem to be ad- 
dressed for the calculation of the Auger theoretical spec- 
trum is the construction of a continuum electronic wave- 
function that accurately takes into account correlation 
effects among electrons bound in the nanoclusters and 
between the escaping and bound electrons. To build this 
wavefunction we need to extend a method [13, [H| already 
successfully applied to the calculation of inner shell spec- 
tra in small molecules, such as CO [29| and C2H2 [30ll3l|. 
However, the existence of several interacting decay paths, 
quantum effects due to complex many-body interactions 
and localization/delocalization effects, make such exten- 
sion very difficult if one wants to keep the computational 
cost comparable to that of molecules. 

To compute the scattering wavefunctions E >, in- 
cluding the appropriate incoming wave boundary condi- 
tions, we use an approach based on Fano's multichannel 
theory 32], where the solution to the scattering process 
is represented by a linear combination of a metastable 
state 1$ > with states embedded in the non- interacting 
continuum many- fold {xp~ e }, 
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with r and M« [ep, E) defined as follows 



r = ^r^27T^|Af^(e /3 ,i?)| 2 (2) 
Mp(ep, E) = - E\ X p.)\ e p = E r - Ep (3) 
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In equations |T][3]) E r is the energy of the resonant state 
modified by the interaction with the continuum states 
{Xg e } and is the energy of final states. These mul- 
tichannel continuum wavefunctions can be obtained from 
a set of interacting continuum states {xp^p}, coupled by 
the interchannel interaction, by projecting the Hamilto- 
nian H on a finite set of multicentered L 2 functions span- 
ning the Hubert space in the range of the ejected electron 
kinetic energy: 

< X pe p \H-E\ Xae „ > = <Tpe f) ;P\ H(E) - E \^ ea ■ a > 

(4) 

where \ip S p] (3 > represents the antisymmetric tensorial 
product of the escaping electron wavefunction at energy 
ep with the doubly ionized target wavefunction (/3). In 
equation (j4|) 7i is the projected Hamiltonian, which in- 
cludes the projected scattering potential 

V^ a ( S0 ,e a )=< X 0e \ V\Xae a > (5) 

describing the interchannel interaction among bound and 
emitted electrons in the scattering region. 

To find an accurate solution of equation ^ in Si0 2 
nanoclusters, a number of numerical and theoretical tech- 
niques have been implemented. The first problem is the 
choice of the basis set for the bound states and the scat- 
tering wavefunction. Due to the large number of atoms 
and open channels in comparison to molecules, the size of 
the Hilbert space where the wavefunctions are expanded 
increases rapidly for nanoclusters. An increasing num- 
ber of Gaussians may lead to convergence to non-physical 
states due to basis set superposition errors and linear de- 
pendence problems [33[ , particularly in the calculation of 
metastable states. In order to keep the number of excited 
determinants in the doubly ionized target low, we split 
the problem in two steps: we first calculate the bound 
orbitals in a restricted Gaussian basis set; we then intro- 
duce an auxiliary basis set of modified Gaussian functions 
to compute the matrix elements of the scattering poten- 
tial (equation ([5])). The Hartree-Fock procedure takes 
the variationally optimized atomic wavefunctions of dif- 
ferent components of the nanocluster as a starting point 
for the calculation of the nanocluster wavefunctions, en- 
suring the convergence of the multiconfiguration proce- 
dure [3J| to the appropriate ionic states within the space 
spanned by a minimal basis set. An adequate representa- 
tion of the continuum orbital, comparable to the solution 
of the multiconfiguration procedure in bound orbitals, 
forces us to enlarge the basis set in scattering wavefunc- 
tion calculations. We emphasize that the scattering basis 
set is built with a diffuse Gaussian basis set multiplied 
by symmetry adapted Hermite polynomials to reproduce 
the oscillating behaviour of the continuum wavefunctions 
inside the scattering region. This enlarged basis set is 
used both in the expansion of the continuum orbital and 
in the projection of the interchannel coupling inside the 
Hilbert space spanned by the interacting decay channels. 
Continuum orbitals have to be orthogonal to bound or- 
bitals. To satisfy this constraint one can use two differ- 



ent methodologies. One is direct orthogonalization, the 
other uses projection of the kinetic operator onto the or- 
thogonal complement of the bound states. The latter 
procedure defines an effective potential, which has the 
same structure of the Phillips-Kleinman potential [35[ 
used in condensed matter calculations for building va- 
lence orbitals orthogonal to the core states. Surprisingly, 
we found that the Phillips-Kleinman technique leads to 
numerical instabilities, with a basis set dependence, in 
the calculation of the Si02 nanoclusters continuum den- 
sity of states. It turned out that such instability is found 
whenever the eigenvalues of the metrics of the chosen ba- 
sis set are smaller than 10~ 4 . Attention has to be paid, 
when choosing the basis set, to describe appropriately 
the scattering potential in the energy region where elec- 
trons are emitted and to avoid a small overlap between 
the basis elements. After taking into account such a cri- 
terion, calculations of the continuum density of states for 
Si02 nanoclusters gave the same results in terms of either 
approach used to satisfy the orthogonality constraint. 

Due to the enormous number of decay channels coupled 
with the final states, the Auger probability in nanoclus- 
ter calculations is spread all over the bands with a com- 
plex distribution. Accordingly, the computational cost 
to evaluate the interchannel potential in equation (JSJ) , re- 
sponsible for the interaction between the open channels 
and for the probability spread among the non interact- 
ing channels, may be so high that it could not be directly 
performed as previously done for molecules [36| . Unfor- 
tunately, interchannel interaction cannot be neglected, 
because its absence is the main source of errors and of 
large discrepancies with experimental spectra in scatter- 
ing calculations [36[. 

One can argue that the strong ionic bond between Si 
and O leads to a rather large interaction. The electronic 
cloud distortion, with silicon to oxygen charge donation, 
favors holes trapping on the oxygen site. Moreover, small 
multisite interaction among the Si02 units, is responsible 
for the narrow valence bandwidth. As a consequence, 
one expects a very low hole mobility, clamped down by a 
large hole-hole interaction. From these observations, one 
expects that the final states are created by intra-atomic 
transitions with the two holes strongly localized. 

In light of what we have shown above, a theoretical 
approach for calculating atomic-like and band-like Auger 
lineshapes in nanoclusters with strong correlation effects 
is proposed here. In this method both the initial state 
core-hole and the final state valence holes are first lo- 
calized on the central atom and multisite interactions 
enter naturally in the model Hamiltonian as a pertur- 
bation. Therefore, a number of localized states, which 
account for the majority of the decay probability, is se- 
lected among the open channels. These 'tight binding' 
states are further coupled via the interchannel potential 
in such a way that the model Hamiltonian can be di- 
agonalized at a very high level of accuracy in a smaller 
space, corresponding to the most intense open channels. 
This calculation results in a number of non interacting 
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states, which account for the quasi-atomic behaviour of 
the nanocluster. The long range or band-like behaviour 
of the Auger decay is calculated at this stage by switch- 
ing on the weak multisite interaction. The perturbative 
coupling of each non-interacting channel with a number 
of states describing transitions close in energy, allows us 
to choose, among the degenerate states, those present- 
ing a maximum overlap with the selected channel and, 
therefore, to obtain the whole lineshape. 

III. CASE STUDY: SI0 2 

To test our approach, we take the Si2 07Hg nanoclus- 
ter as representative of bulk SiC>2 to calculate the normal 
O K-LL Auger spectrum. Within the Born-Oppenheimer 
approximation, the optimized structure of Si207Hg with 
a central oxygen atom has been obtained by minimiz- 
ing the free energy of the electrons at each nuclei posi- 
tion with a smooth Fermi-Dirac smearing, relaxing the 
atomic bond length until the forces are less then 0.01 
eV/A per atom. The resulting atomic configuration is 
plotted in Fig. (TT]). Stars of neighbors, centered on the 
central oxygen atom, are found at a distance of 3.103 ao 
for silicons and 4.953 ao for oxygens. The present re- 
sults are in good agreement with previous calculations 
on band structure and geometry of similar Si02 clus- 
ters [l?]]- Following our scheme, correlation in the ini- 




FIG. 1: '(Color online)' Sketch of the Si02 nanocluster op- 
timized structure. Oxygens are in red, silicons in light blue, 
hydrogens in green. 

tial state of the ionized target, characterized by a single 
inner vacancy in the Is orbital of the central oxygen, 
and in the final states of the doubly ionized nanocluster, 
characterized by two vacancies in the valence orbitals, 
has been treated by extended configuration interactions 
(CI) expansion using the 3s, 3p, 3d orbitals of silicon and 
3s, 3p, 3d of oxygen. Calculations of the bound states 
have been carried out using a basis set of Hermite Gaus- 
sian functions with (s, p, d)-type character centered on 



TABLE 1: Kinetic energies (Ekin) in eV, referred to the vac- 
uum level, and probabilities (r a ) in arbitrary units for the O 
K-LL Auger localized states of Si207H6 nanocluster according 
to the double holes configu-rations (O K-LL) in the central 
oxygen atom and final total spin S 2 (0=singlet, l=triplet). 
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the nuclei. The basis set has been obtained by a varia- 
tional procedure, optimizing the exponents of the stan- 
dard 6-31G* basis set. In more detail we scan a finite 
range of values above and below the proposed 6-3 1G* 
ones, taking the values which lower the total energy. An 
auxiliary basis set, more dense and enlarged with diffuse 
functions, has been used to represent properly the pro- 
jected potential V 3 '" in equation (5), until variational 
stability of the Auger decay rate with respect to changes 
in the basis set is reached. This criterion allows an ac- 
curate description of the continuum orbitals around the 
Auger decay energies. The scattering calculations high- 
light the fundamental role of d-symmetry functions in 
the expansion of the projected potential. The resulting 
basis set is made up of [9s(= Assi + 3so + 2sh) + 6p(= 
3psi + 2po + lpif) + 3d(= l^Si + 2do)] contracted Her- 
mite Gaussian functions, for a total of 215 contracted 
Gaussians. The exponents a s ,a p ,ad are in the range 
[af = 0.001^-1.0, af = 0.004-=- 1.2, af = 1.0] for the (s, 
p, d)-type Gaussians on silicons, [af = 0.002 -~ 0.5, a° = 
0.07-=- 1.0, af = 1.0] for the (s, p, d)-type Gaussians on 
oxygens, [af = 0.03 -j- 1.0, a p = 1.0] for the (s, p)-type 
Gaussians on hydrogens. We adopted a larger basis set 
on the central oxygen with [4s + Ap + 3d] to properly 
represent hole localization effects. The examined nan- 
ocluster has a ionic character, with the central oxygen 
atom in the Si20 unit found in the Si +149 0~ 120 charge 
configuration. We found that the Si207H6 nanocluster 
correctly reproduces one electron properties, localization 
effects and the charge population of solid Si02- In table 
(fljl we have reported the energies, referred to the vacuum 
level, and the relative decay rates of the most intense fif- 
teen transitions, all with double holes localized on the 
central oxygen atom, obtained after switching on the in- 
terchannel coupling. 

IV. EXPERIMENTAL DATA HANDLING 

Since first-principles methods necessarily involve ap- 
proximations, the most important of which, in scatter- 
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ing calculations, is the approximate form of the electron 
wavefunction in the continuum, it is important to vali- 
date the calculations against experimental data. To this 
end, O K-LL spectra were measured from a 36 nm thick 
SiC>2 film, grown on a Si substrate. The film surface was 
cleaned by annealing at 800° C in a Ultra High Vacuum 
preparation chamber connected to the analysis cham- 
ber (base pressure 2xl0 -10 mba). The analysis cham- 
ber is a PHI545 instrument equipped with double-pass 
Cylindrical Mirror Analyzer, coaxial electron gun, non- 
monochromatic Mg Ka X-ray source and He discharge 
lamp. The X-Ray (tn/=1253.6 eV) excited O K-LL spec- 
trum was recorded at constant analyzer energy resolu- 
tion, AE=0.6 eV, as measured on the Pd Fermi edge of 
a Hel (h^=21.2 eV) excited valence band photo emission 
spectrum. A Reflection Electron Energy Loss Spectrum, 
excited from the same sample by 500 eV electrons, i.e. 
near the kinetic energy of O K-LL electrons, was acquired 
within the same instrument and at the same analyzer res- 
olution of the Auger spectrum. To keep electron beam 
induced damage to a minimum, the current density was 
lower than 40 A/m 2 and the acquisition time was lower 
than 6 minutes. The energy scale, was calibrated on clean 
Au and Cu, following recommended procedures [38| . Ki- 
netic energies are referred to the vacuum level. 

'As acquired' spectra were corrected for the energy de- 
pendence of the analyzer transmission function. The lat- 
ter had been previously determined on clean Au, Ag and 
Cu as explained in [391 ]. A constant background, defined 
by the region immediately to the right (520-530 eV) of 
the main peak in the O K-LL spectrum, was subtracted, 
and intensity was normalized to unit height of the main 
peak, resulting in the spectrum (dotted line) shown in 
Fig. ©. Also plotted in Fig. © is the REEL spec- 
trum (dashed line) which includes the so called elastic 
or zero loss peak together with the associated inelastic 
losses. The spectrum is slightly shifted in energy, so that 
the zero loss peak is aligned to the main O K-LL peak. 
Its intensity is normalized to unit height of the zero loss 
peak. A vertical expansion of the loss features is shown 
in the inset. By deconvoluting the REEL spectrum from 
the measured O K-LL spectrum, the O K-LL line-shape 
given by the continuous line in Fig. ([2]) was finally ob- 
tained. 



V. LINESHAPE ANALYSIS 

This spectrum is compared to the ab-initio calculated 
O K-LL spectrum in Fig. ([3]). The computed spectrum is 
generated by convoluting the theoretical transition rates 
with a Voigt profile to take into account both the band- 
width of the X-ray source and the finite resolving power 
of the electron spectrometer. This results in a convolu- 
tion of a Gaussian with S g = 1 eV, according to our in- 
strument resolution, and a Lorentzian with S; = 0.1 eV, 
which is the computed total Auger decay rate. Both spec- 
tra are normalized to a common height of the main peak, 
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FIG. 2: Plot of the experimental data. The X-Ray 
(h^=1253.6 eV) excited O K-LL spectrum (continuous line) 
has been obtained by deconvolution of a REEL spectrum 
(dashed line) from a spectrum obtained subtracting a con- 
stant background, defined by the region immediately to the 
right (520-530 eV) of the main peak (dotted line). 



while no energy shift was required to obtain the kind of 
agreement shown in the figure. We observe that our cal- 
culation correctly reproduces transition energies. These 
are distributed over three different regions. From right 
to left, features between [494 — 510] eV are due to 2p — 2p 
hole configurations in the final states (K-L23L23 tran- 
sitions); the most intense peaks are singlet transitions 
around [501 — 503] eV (maximum at 502.76 eV), close to 
the experimental peak (502.58 eV). Between [475 - 492] 
eV a second group of features is found, where the hole 
configuration in the final states is 2s — 2p (K-L1L23 tran- 
sitions); the most intense peaks are singlet transitions 
around [480 - 482] eV (maximum at 481.1 eV), close to 
the experimental peak (481.38 eV). Finally a third group 
of features is observed between [460 — 470] eV, where the 
hole configuration in the final states is 2s — 2s (K-L1L1 
transitions); the most intense peaks are singlet transi- 
tions around [463 — 465] eV (maximum at 463.86 eV), 
lower than the experimental peak (465.25 eV). In general 
the most intense peaks are singlet transitions with both 
holes on the central oxygen atom, whereas the shoulders 
are triplet transitions, unfavored by selection rules. 

Besides transition energies, also calculated values of 
the decay rates are satisfactorily reproduced in Fig. ([3]), 
meaning that the adopted scattering basis set is appropri- 
ate to represent the projected potential within the scat- 
tering region for any channel. However, the ss and sp con- 
tributions are slightly overestimated. This is due to lack 
of electronic correlation or to dynamic screening effects 
[Tij in the CI calculation. The size of the excited determi- 
nants space in the CI procedure for both the intermediate 
core-hole and double ionized final states has been chosen 
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large enough to reproduce well transition energies and 
decay probabilities, but a larger active space should be, 
in principle, adopted to take into account hole dereal- 
ization on all valence bands and electronic excitations. 
Furthermore, the localized (not atomic) Auger matrix 
elements are kept fixed during this procedure and their 
square modulus enters as a multiplicative factor in the 
lineshape analysis. In principle, they should be updated 
with a new electronic continuum wavefunction derived by 
a multichannel procedure, taking into account the cor- 
relation between the double ion and the new escaping 
electron wavefunction. Larger CI calculations, though 
computationally very expensive, would correct the line- 
shape for the overestimation of the s to p contribution. 
The same kind of agreement seen in Fig. (|3|) is obtained 
by comparing our calculated O K-LL Auger spectrum to 
the measured spectrum by van Riessen et al. [8| and Ra- 
maker et al. [l(|. Without loss of accuracy the present 
method thus allows us to calculate the Auger lineshape in 
nanostructures and solids at a computational cost com- 
parable to that typical for atoms and molecules. 
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FIG. 3: O K-LL Auger spectrum in Si2 07H6. Comparison be- 
tween the Quantum Mechanical (QM) theoretical data (dot- 
ted line) and the experimental data after deconvolution of a 
REEL spectrum (continuous line). 



probability distribution as a source of electrons which 
undergo inelastic processes. Within our model the ab- 
initio spectrum could include the intrinsic energy loss, 
due to shake processes. As stated above, however, shake 
processes were neglected in the present paper. On the 
other hand, extrinsic energy losses, due to the interaction 
with collective degrees of freedom of the valence charge 
(valence plasmons [4(|) and to multiple inelastic scatter- 
ing, should be dealt with after computing the spectrum 
and using statistical methods. In this way inelastic back- 
ground would be superimposed to the theoretical spec- 
trum instead of removing it, a procedure not free from 
uncertainty. Possible effects of background subtraction 
on the spectrum final shape are discussed, for example, 
in 

The simplest approach to account for the extrinsic loss 
is the Monte Carlo method [42|. Details of the Monte 
Carlo scheme describing electron energy losses in sili- 
con dioxide were given in previous works (43l . IlH . |45| . 
In the following we will briefly summarize the numerical 
approach used in this paper. The theoretical Auger spec- 
trum previously calculated is assumed to be the initial en- 
ergy distribution of the escaping electrons. The stochas- 
tic process for multiple scattering follows a Poisson-type 
law. The step-length As is given by As = — A ln(/ii) , 
where \i\ is a random number uniformly distributed in 
the range [0, 1]. A is the electron mean free path: 



X(E) 



1 



N[a el {E) + a mel (E)} 



(6) 



where N is the number of Si02 molecules per unit vol- 
ume, <J e i{E) is total elastic scattering cross section, and 
Cmei(E) is the total inelastic scattering cross section at 
the impinging electron kinetic energy E. The calculation 
of differential and total elastic scattering cross sections 
was performed by the Relativistic Partial Wave Expan- 
sion method [43|, |4(| [47|, whereas differential and total 
inelastic scattering cross sections were computed using 
the Ritchie theory (48|. Taking Buechner experimental 
optical data [4!| for the evaluation of the long-wavelenght 
limit of the dielectric function e(u>), one obtains the dif- 
ferential inelastic scattering cross section as: 



da ine i(E,uj) 
dui 



2irh 2 NE 



Im 



" -1 " 











(7) 



VI. MONTE CARLO SCHEME: QMMC VS. 
MEASURED SPECTRUM 



As stated above, the comparison between computed 
and experimental spectra requires to account for changes 
caused to the original electron energy distribution by en- 
ergy losses suffered by the Auger electron on its way out 
of the solid. The typical procedure in this respect is the 
one described above. Alternatively however, one could 
simulate the effects of inelastic losses on the 'original' 
distribution. The basic idea is to use the ab-initio Auger 



where u> is the energy loss and, according to Ashley [5C 
the function S is given by 



S{x) = (l-x)ln- 



7 



,3/2 



33 
32 : 



(8) 



Before each collision, a random number /12 uniformly 
distributed in the range [0, 1] is generated and com- 
pared with the probability of inelastic scattering qi ne i — 
Cinei I (&inel + &el) ■ When fi2 is smaller or equal to the 
probability of inelastic scattering, then the collision is in- 
elastic; otherwise, it is elastic. If the collision is elastic, 
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9, the polar scattering angle, is selected so that the inte- 
grated scattering probability in the range [0, 9] is equal 
to the random number /13 uniformly distributed in the 
range [0,1]: 

Ms = — / ^2tt smddti, (9) 

Pel Jo 



VII. CLUSTER SIZE EFFECTS 

Since we are dealing with nanoclusters, it makes sense 
to assess the dependence of correlation effects on the clus- 
ter size. It is well known that an increasing number of 
atoms in the cluster causes the one-electron DOS to ap- 
proach the band DOS of the solid. Increasing the number 



where ft is the solid angle of scattering. If the collision is 
inelastic, the energy loss W of an incident electron with 
kinetic energy E is computed via a random number /X4, 
uniformly distributed in the range [0,1]: 



1 



&inel JO 



ir 



(III! 



■du> 



(10) 



Auger electron generation is simulated assuming a con- 
stant depth distribution whose thickness, according to 
Q, was set to 40 A. To obtain good accuracy in the re- 
sults the number of random walks is 10 8 . A plot of the 
QMMC calculation compared to the experimental data 
before deconvolution of energy losses is given in Fig. 
The original theoretical spectrum is also shown for refer- 
ence. One can see that QMMC enhances and broadens 
the Auger probability increasingly upon decreasing the 
electron kinetic energy. The large broadening of the K- 
LiL2,3 peak after Monte Carlo treatment is due to the 
main plasmon of Si02 (see inset in Fig. whose dis- 
tance from the zero loss peak is the same as the distance 
between the K-L23L23 and K-L1L23 features in the Auger 
spectrum. 




460 480 500 
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FIG. 4: O K-LL Auger spectrum in Si2C>7H6. Comparison 
between the Quantum Mechanical theoretical data (dotted 
line), the Monte Carlo results (dashed line) and the original 
experimental data (continuous line) . 
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FIG. 5: Comparison between the O K-LL spectrum in Si20He 
(dotted line), Si207Hg (continuous line) and experimental 
data (dashed line). 

of atoms to simulate a bulk situation enlarges the number 
of available transitions to the final states to a point where 
the computational cost is enormous, while no appreciable 
change in the spectroscopic observables is obtained. In 
fact, the quasi-atomic nature of Auger emission found in 
our ab-initio calculations forces the hole-hole DOS to be 
centered on the atom where the initial core-hole has been 
created. We have investigated the role of cluster size by 
carrying out calculations on a smaller cluster of atoms, 
Si20He, to see how it compares with the bigger clus- 
ter. This nanocluster represents the minimal cluster size 
which still describes the chemical environment found in 
solid Si02- In Fig. ([5]), we compare the calculated O K- 
LL spectrum of Si2 0H 6 (dotted line) to that of Si2 07H 6 
(continuous line). The experimental spectrum is also su- 
perimposed (dashed line). All spectra are normalized to a 
common height of the main peak, while no energy shift is 
performed. The two computed spectra were obtained by 
convolution with the same Voigt profile. We see that the 
main features of the spectrum are already consistently 
reproduced in the smaller cluster, though with a shift to- 
wards high energies. This 'chemical shift' decreases from 
2.52 to 1.45 eV on moving from the low energy to the 
high energy peak. This is a predictable consequence of 
the fact that Si dangling bonds are saturated with differ- 
ent elements for the two clusters (H for the smaller cluster 
and O for the bigger one) with a non-local dependence 
which is larger for lower angular momenta. Moreover, we 
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observe that the Si2 07Hg lincshape is less resolved than 
the corresponding Si20H 6 one, as expected due to the 
increased number of decay channels. Though the two 
computed spectra do not differ much in the number of 
features and in relative decay rates, the overall agree- 
ment with the experimental spectrum slightly favors the 
bigger cluster. This observation confirms that our cluster 
approach is well-suited for the calculation of Auger line- 
shapes from bulk materials and that the f 5 atom cluster 
considered in the present paper is near to the optimal 
size. In conclusion, we believe that enlarging further the 
cluster would not change substantially the Auger spec- 
tral lincshape as expected on the basis of the local nature 
of AES. Inclusion, in the computed spectrum, of shake- 
up and shake-off transitions and of interactions between 
crystal phonons and escaping electron would result in a 
broadening of the lineshape, correcting for the difference 



between the computed and experimental linewidth. 



VIII. CONCLUSIONS 

We have demonstrated the potentiality of our first 
principle model for calculating and predicting the line- 
shapes of X-Ray excited Auger Electron Spectra from sil- 
icon dioxide nanoclusters. Furthermore, a Monte Carlo 
analysis has been applied to the Quantum Mechanical 
computation of spectra to mimic the electron energy loss 
in condensed matter. Comparison to our experimental 
data reveals a satisfactory agreement, in the energy po- 
sition and relative intensities of peaks, and in the back- 
ground contribution over the entire investigated energy 
range. 
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